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چکیده امروزه بسیاری از مناب عآب‌های سطحی و زیرزمینی در اثر ورود پساب‌های ناشی از مواد مختاف در معر ضآلودگی قرار می‌گیرند. با توجه به اینکه اکثر 
روش‌ها جهت بازیابی تابع شد تآلودگی برا ی آب‌های زیرزمینی انجا مگرفته است, ارائه روشی جهت بازیابی تابع شدت آلاینده در رودخانه‌ها با شرایط جریان 
پیچیده‌تر می‌تواند مورد توجه قرا رگیرد. در این پژوهش هدف محاسبه تابع شدت منبعآلاینده بوده که بر اساس تحقیقات صورت گرفته تاکنون این مسئله با 
استفاده از رامحل عددی طرح حافظ گروه مورد مطالعه قرار نگرفته است. طرح حافظ گروه روشی دفیق برای حل مسائل بدخیم است که در این نوشتار با استفاده 
از این روش, معادله یک‌بعدی جابه‌جایی-پراکندگی با ضرایب متغیر حل شده است. اساس روش حل معکوس ارائه شده به این صورت است که با حل سیستم 
های دینامیکی در گام‌های زمان‌های مفیء یک معادله کلی حاصل می‌شود که قادر به حل معادلات دیفرانسیل معمولی حواهد بود. پاسخ‌های ب رگشتی این معادله 
را به همگرایی لازم می‌رساند و از واگر! شدن داده‌ها جلوگیری می‌نماید. در اين پژوهش سه مثال جهت نشان دادن دقت طرح ستقیم و معکوس حافظ گروه 
ارائه شده است. ابتدا با استفاده از حل مستقیم مقدار غاظت آلاینده در رودخانه با ضرایب متغیر محاسبه شده است و به‌منظور صحت‌سنجی حل مستقیم از حل 
عددی تمام ضمنی تفاضل محدود استفاده شده است. مقایسه نتایج این دو روش با شاحص‌ها یآماری, انطباق این دو مدل در روش مستقیم را نشان می‌دهد. در 
گام بعد با دو ال مختلف با استفاده از حل طرح معکوس حافظ گروه به شبیه‌سازی تابع شدت آلاینده در بلادست رودخانه پردانعته می‌شود و پس از شبیه‌سازی 
و صحت‌سنجی مدل معکوس با حل مستقیم, ارزیابی این روش با استفاده از شاحص‌ها یآماری انجام می‌شود که نتایج دقت خوب این روش را نشان می‌دهد. 
واژه‌های کلیدی حل عددی. روش معکوس. طرح معکوس حافظ گروه معادله جابه‌جایی- پراکندگی, روش عددی خطوط. 
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مقدمه 


امروزه آلودگی آب‌های سطحی یکی از بزرگترین معضلات 
زندگی انسان و از مهم‌ترین مسائل زیست‌محیطی است که 
برخلاف اهمیت آن در معرض انواع آلودگی‌های شیمیایی و 
بیولوژیک قرار دارد. باوجود قوانین و مقررات برای جلوگیری از 
تخلیه غیرمجاز پساب‌های ناشی از صنایع و کشاورزی. بازهم 
جلوگیری از اين اتفاق امکان‌پذیر نیست. ازاین‌رو شناخت مکان 
و زمان تخلیه آلودگی در کمترین زمان و کمترین هزینه موردتوجه 
محققان این بخش قرار گرفته است. روش‌های مسئله معکوس 
(صعاام2 6وت2۲7۵]) رامحل بسیاری از مشکلات در شناسایی 
ی و منبع آلاینده می‌باشد. با توجه به اهمیت حفاظت از منابع 
آب در دهه‌های اخیر, شناسایی منابع آلاینده موردتوجه بسیاری 
ز پژوهش گران قرار گرفته است. 

می‌توان گفت مهمترین پارامتر انتقال آلودگی در رود خانه 
نتقال آلودگی از یک مکان به مکان دیگر است که توسط جریان 


ب‌ رخ 7 هد. این فرایند اند ال جابه جایی (ظ۸0۷6۵۱۲0) 
من‌تاشته 1 از نکر فرانتتاعال فد رها نلفر ایند بزا کند کر 
(1(150675100) می‌باشد که توسط فیک (۳0) معرفی شده است 
11 در بحث انتقال آلودگی در رودخانه عموماً فرایند جابه‌جایی 
مطرح اسست اما فرایند پراکندگی نیز اثرگذار می‌باشد [2]. بر 
اساس قانون اول فیک. معادله دیفرانسیل حاکم بر انتقال و انتشار 
آلاینده در رودخانه همان معادله جابه جایی-پراکندگی-واکنش 
(ممتاعدم-۸0۷6۵0۲-1۲15069102) است که از نوع معادلات 
دیفرانسیل سهموی غیر همگن, خطی و مرتبه دوم می‌باشد و با 
ها رس و ی کم فا ای سم زان گرد 
شدت منبع آلاینده را با توجه به معلوم بودن مکان و زمان توزیع 
غلظت آلودگی محاسبه نمود. 

به‌طورکلی حل معکوس معادله انتقال آلودگی در رودخانه‌ها 
و آب‌های زیرزمینی توسط سه روش کلی شبیه‌سازیبهینه‌سازی 
(0مطاعصه ممتاهعتهتامه صمتاداصن8). احتمالاتی و ز مین‌آمار ی 
(طمعمتووه لمع‌تادناداد0ع0)) و ریاضی (0مط)عصه 21عت0امصصعط۷]۵) 
مورد بررسی قرار گرفته است [3]. 
در روش اول از ترکیب یک الگوریتم بهینه‌سازی با یک روش 
عددی جهت حل معکوس معادله انتقال آلودگی استفاده می‌شود 
که به طبع نیاز به کامپیوترهایی با پردازنده قوی می‌باشد که نتیجه 
آن افزايش هزینه و زمان محاسبات می‌باشد و از مزیت‌های این 


نشريهةُ مهندسی عمران فردوسی 


حل معکوس (در زمان) معادله انتقال آلودگی در رودنحانه... 


روش می‌توان به‌تناسب این روش با محیط‌هایی با مقیاس بزرگ 
آشاره کردر هر هم انا دز سال ۳۰۱۱۵ تاردو و عمکاران ]نا 
استفاده از روش تنظیم حداقل مربعات ( ععتقناو؟ 15 
6 و همچنین استفاده از اطلاعات دریافت شده از 
حسگرهای زیستی, به شبیه‌سازی زمان و منبع آلاینده در محیط 
یک‌بعدی پرداختند. در سال ۲۰۱۷ نیز ژانگ و ژین [5] با ترکیب 
روش‌های الگوریتم ژنتیک پایه و حل تحلیلی معادله انتقال و 
انتشار تابع شدت زمانی آلاینده را بازسازی کردند. مهینتاکومار و 
تیار سا ۱۱ وه شک ۱۱۱7۱۰۱۳۳ اشفا 
از این رون به شیه‌سازی مادله اتقال اتسار پرداغشند 161 

در روش دوم از کاربرد روش‌های احتمالاتی و زمین‌آماری 
در شبیه‌سازی مکان و زمان غلظت آلاینده استفاده می‌شود. از 
مزیت‌های این روش می‌توان به کاهش حجم محاسبات و از 
معایب آن می‌توان به محدودیت در شناسایی تعداد منابع آلاینده 
اشاره کرد. در همین راستا در حل معکوس معادله انتقال آلودگی 
در رودخانه با استفاده از روش احتمال برگشتی ( 501270 
۵مطع نآنطدط۳:0) تحقیقاتی توسط پژوهشگرانی انجام شد. 
در سال ۲۰۱۹ قانع و همکاران[8] با استفاده از این روش محل و 
زمان انتشار منبع آلاینده را محاسبه کردند و نشان دادند که این 
روش دقت محاسباتی مناسبی دارد. همچنین این پژوهشگران در 
پژوهشی دیگر دقت خوب این روش در شناسایی زمان انتشار و 
محل منبع آلاینده در رودخانه را نشان دادند. در پژوهشی ۳ 
چنگ و جیا [9] در سال ۲۰۱۰ با استفاده از روش احتمال برگشتی 
در پهنه‌های آب‌های سطحی و با شرایط جریان غیریکنواخت و 
ماندگار و در حالت دوبعدی موفق به تشخیص مکان و زمان 
رهاسازی ماده آلاینده در رودخانه شدند و با صحت‌سنجی این 
روش با داده‌های واقعی دقت خوب این روش را نشان دادند. 
نیوپار و ویلسون [10] در سال ۲۰۰۵ با روش احتمالات برگشتی. 
کاملا وهنکایان |۱1 سال ۱۱۵ نوی یی بر وهای 
زمینآماری و در سال ۲۰۱۸ وانگ و همکاران [12] با استفاده از 
روش احتمال برگشتی و مدل رگرسیون خطی ( 10627 
۸۵ ۲68۲6۹510۴) در همین راستا تحقیقاتی انجام دادند. 

در روش سوم با استفاده از چارچوب‌های ریاضی شبیه‌سازی 
انجام می‌گیرد. روش‌های ریاضی به دلیل کاهش زمان محاسبات؛ 
کاهش تعداد دفعات اجرا مدل و حذف کامل تکرارها مزیت 


بیشتری نسبت به روش شبیه‌سازی-بهینه‌سازی دارد و از معایب 


سال سی و پنجم» شمارة چهار ۱۶۰۱ 


امیر محمد سعادت- مهدی مظاهری - جمال محمدولی سامانی 


این روش می‌توان به پیچیدگی حل مسئله نسبت به سایر روش‌ها 
اشاره کرد. در حل معادله انتقال و انتشار با استفاده از روش 
ریاضی می‌توان به حل معکوس معادله انتقال آلودگی در رودخانه 
با ضرایب متغیر با استفاده از روش شبه معکوس‌پذیری اشاره 
کرد. لوشابی در پژوهشی به حل معکوس معادله یک‌بعدی انتقال 
آلودگی با ضرایب متغیر پرداخته است که با استفاده از حل شبه 
معکوس پذیری معادله جابه‌جایی-پراکندگی و یک مدل عددی؛ 
داده‌های سری زمانی غاظت آلاینده یفده کر وودخانه 


متسه ق باسازی شنم اسها |۲9 خرسان ۳۲۰ تبرت 1241] 
نیز در پووهشی به حل معادله تال آلودگی با استخراج راهحل 
تحلیلی توسط روش شبه معکوس‌پذیری پرداخت که با افزودن 
ترم پایداری شبیه‌سازی معادله جابه‌جایی-پراکندگی را انجام داد. 
نتایج نشان داد که این روش در الگوهای ساده بسیار دقیق عمل 
یکلا ور هن الکو ها میم ان قطا اف آیش.می بانن: با انتفاده 
از اين روش پژوهش‌هایی توسط اسکاگاس و کابالا [15] در سال 
۵ وه زانی و فانک [16]کد سال ۲۰۱۱ به شبه‌شارعن سادله 
انتقال آلاینده پرداختند. از روش‌های عددی و تحلیلی مختلف 
دیگری همچون روش تنظیم تیخنوف. روش تابع گرین تبدیل 
لاپلاس و غیره جهت حل معکوس معادله انتقال آلودگی در 
زونه بش ایط معلت: انتفاوه شله استا که دز ین ,راما 
ای کار ان نی تیال ۲۵۸ نا اناوت او ره فش توف 
(60000 صمتام2تتماناعع: «محمط۲11) . مظاهری و همکاران در 
سال ۲۰۱۵ با استفاده از معادله انتگرالی حاصل از اعمال روش 
تابع گرین به شبیه‌سازی مکان و زمان آلاینده در رودخانه 
پرداختند [3,17] 

با توجه به اينکه در اين پژوهش برای حل معکوس معادله 
انتقال آلودگی از طرح معکوس حافظ گروه ( ونام6) ۳0۱/274 
معط ع۳:۵۹۵۲۳۷[9) استفاده شده است. در اینجا به‌صورت 
مختصر به تحقیقات مختلف که توسط طرح حافظ گروه (0:0۲ 
م56 ۳۳۵5۵۲۷12۵) صورت گر فته پرداخته می‌شود. 

لئو[18] در سال ۲۰۰۶ با استفاده از طرح حافظ گروه به حل 
معکویتن وله ال بجر ریت داسی د اتن مضل با اسکفاده از 
شرایط مرزی دیریکله (1017100166) در یک دامنه فضایی مشخص. 
ابتدا با استفاده از روش نیمه گسسته‌سازی به جداسازی معادله 
انتقال حرارت پرداخت و پس از تبدیل معادله دیفرانسیل جزئی 
به معادله دیفرانسیل معمولی. با استفاده از طرح معکوس حافظ 
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گروه و با حل سیستم‌های دینامیکی به یک معادله کلی دست‌یافت 
که با استفاده از آن به حل معادله دیفرانسیل به‌دست‌آمده در حالت 
یک‌بعدی و دوبعدی پرداخت و نتایج روش طرح معکوس حافظ 
گروه را برای شرایط مختلف سنجید و رنج پایداری را برای این 
معادله تعیین کرد. در سال ۲۰۱۰ لیو و چانگ[19] با فرض ثابت 
بودن ضریب پراکندگی در نواحی متناظر و با استفاده از روش 
عددی خطوط و طرح معکوس حافظ گروه به‌منظور شناسایی 
ویژگی‌های منبع آلودگی شامل محل منبع آلاینده. بازیابی 
تاریخچه انتشار و توزیع فضایی غلظت آلاینده در فضای یک- 
بعدی به حل معکوس معادله انتقال آلودگی در آب‌های زیرزمینی 
پرداختند. همچنین دو روش طرح معکوس حافظ گروه 
و۳[ (صمتافنمظ صحعظ هسامح رها -ممتطه‌ته]) را 
باهم مقایسه کردند که نشان دادند طرح معکوس حافظ گروه در 
حل عددی معادله انتقال آلودگی و گرما نتایج بسیار دقیقتری به 
دست می‌دهد. در پژوهشی دیگر که در سال ۲۰۱۶ چانگ و 
کو[20] انجام دادند با استفاده از طرح معکوس حافظ گروه به 
حل معادله غیرخطی کلاین -گوردون پرداختند. در این پژوهش 
آنها معادله بد وضع و دوبعدی کلاین-گوردون ر.» پس از 
جداسازی با شرایط مرزی معین و تبدیل آن به معادله دیفرانسیل 
معمولی با استفاده از طرح معکوس حافظ گروه حل کردند و 
سپس پایداری این روش را موردبررسی قراردادند و نتایج بسیار 
مطلوبی از روش طرح حافظ گروه به دست آوردند. در سال 
اش تن ماک 211 مقادله ای با که از 
معادلات دیفرانسیل جزئی بدخیم است با استفاده از روش نیمه 
گسسته‌سازی به معادله دیفرانسیل معمولی تبدیل کردند و سپس 
با استفاده از طرح حافظ گروه به بررسی و حل معادله پرداختند. 
نتایج کار موفق بودن ترکیب روش خطوط و طرح حافظ گروه 
برای به دست آوردن جواب معادله لاپلاس بدخیم را نشان داد. 
در پژوهشی دیگر گاوو و همکاران [22] در سال ۲۰۱۹ قدرت 
طرح حافظ گروه در شبیه‌سازی مدل کلاین -گوردون که از مسائل 
بد وضع می‌باشد را نشان دادند. آنها به کمک روش عددی خطوط 
و تبدیل معادله دیفرانسیل جزئی به معادله دیفرانسیل معمولی و 
شبیه‌سازی آن با استفاده از طرح حافظ گروه کارآمد بودن اين 
روش در حل معادله کلاین-گوردون در حالت یک‌بعدی و 


دوبعدی را نشان دادند. 
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چانگ و لیو[23] در سال ۲۰۱۶ با استفاده از طرح معکوس 
حافظ گروه در محیط ناهمگن, به حل معکوس معادله چندبعدی 
و غیرخحطی موج (مصواهام0۲ ۸/۵۷۵ 0201607۷2۲0) که جزو معادلات 
دیفرانسیل بدخیم است پرداختند. آنها پس از اعمال شرایط مرزی 
و جداسازی به روش نیمه گسسته‌سازی و تبدیل معادله موج به 
معادله دیفرانسیل معمولی, با استفاده از طرح معکوس حافظ گروه 
به حل آن پرداختند. در ادامه با حل مثال‌هایی پایداری و نتایج 
این روش بر روی معادلات چندبعدی و بدخیم را بررسی کردند. 
بر این اساس با مثال‌های عددی نشان دادند که طرح معکوس 
حافظ گروه برای سیستم‌های غیرخطی چندبعدی و حتی برای آن 
دسته از معادلات که به‌ شدت بدخیم هستند قابل استفاده است و 
همچنین نشان دادن که طرح معکوس حافظ گروه می‌تواند نتایج 
صحیح و دقبقی را تولید کند و پیاده‌سازی این روش در شرایط 
زمانی معین به‌خوبی قابلاجرا است. در سال ۲۰۱۸ شن [24] با 
استفاده از طرح حافظ گروه به حل معکوس معادله غیرخحطی و 
ناهمگن انتقال حرارت پرداخت که جزو معادلات دیفرانسیل 
بدخیم می‌باشد. شن برای تبدیل این معادله به معادلات دیفرانسیل 
معمولی از جداسازی به روش نیمه گسسته‌سازی استفاده کرد و 
با استفاده از طرح حافظ گروه به حل آن پرداخت و همچنین با 
مغال‌های. تحت قترایط متختلفت: بانذاری: :و جواب‌ها زا پزرسی 
کرد و نتیجه گرفت که طرح حافظ گروه بسیار دقیق, پایدار مژثر 
و حساس به داده‌های نهایی هستند و همچنین سرعت‌های 
همگرایی سریعی دارند. روش‌هایی نیز برای تخمین ضریب 
اکن کون تزا تشه اش کر هش اشامن توان به 
پژوهشی که نوری و همکاران[25] در سال ۲۰۱۷ انجام دادند 
اشاره کرد. در این پژوهش با استفاده از مدل محاسبه دانه‌ای که 
بر مبنای روش عددی می‌باشد. به تخمین ضریب پراکندگی طولی 
در رودخانه پرداختند. در پژوهشی دیگر در سال ۲۰۱۹ رمضانی 
و همکاران[26] نیز با استفاده از راه حلی بر مبنای تفاضل محدود 
روشی برای محاسبه ضریب پراکندگی ارائه کردند. 

با توجه به پژوهش‌هایی که تاکنون انجام گرفته. حل معکوس 
معادله انتقال آلودگی در رودخانه با ضرایب متغیر به‌صورت 
محدود صورت گرفته است و طبق بررسی‌های صورت گرفته 
هیچ مطالعه‌ای در خصوص روش شناسایی تابع شدت منبع 
آلاینده در آب‌های سطحی با استفاده از طرح معکوس حافظ گروه 
انجام نشده است. پس در این تحقیق سعی بر آن شده تا با اعمال 
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حل معکوس (در زمان) معادله اتقال آلودگی در رودنحانه... 


طرح معکوس حافظ گروه برای حل معادله جابه‌جایی-پراکند گی 
در حالت یک‌بعدی و ضرایب متغیر در رودخانه برای شناسایی 
تابع شدت منبع آلاینده از مزایای این روش استفاده کافی شود. 
از مزیت‌های این روش می‌توان به عدم استفاده از ترم پایداری 
در حل معکوس معادله انتقال آلودگی اشاره کرد و مشکلاتی 
همچون انتخاب و بهینه‌سازی ترم پایداری را ندارد. مزیت دیگر 
این روش برخلاف ساير روش‌های عددی» عدم نیاز به 
تقسیم‌بندی رودخانه به بازه‌های زیاد می‌باشد که با توجه به‌وقتگیر 
بودن و هزینه زیاد داده‌برداری این نکته اهمیت زیادی دارد. در 
واقع همگرا شدن با تعداد بازه‌های کم یا گام زمانی یکی از 
محاسن این روش است که در این پژوهش نیز این امر اتخاذ شده 
است. سرعت و دقت بالای این مدل از دیگر مزیت‌های مهم 


است که می‌توان به آن اشاره کرد. 


مبانی تثوری و روش‌ها 


در این پژوهش به کمک طرح حافظ گروه که روشی بر پایه 
گروه‌های لی می‌باشد به حل مستقیم و معکوس معادله انتقال 
آلودگی در رودخانه پرداخته می‌شود. در این بخش ابتدا طرح 
حافظ گروه که جهت حل مستقیم معادله انتقال آلودگی کاربرد 
دارد معرفی گردیده و در ادامه به معرفی طرح معکوس حافظ 
گروه که برای حل معکوس این معادله کاربرد دارد پرداخته می - 
شود. در انتها معادله جابه‌جایی-پراکندگی به‌صورت مستقیم و 
معکوس با استفاده از این زوش حل‌شده است. در این پژوهش 
جهت حل مسئله ابتدا با روش عددی مستقیم معادله انتقال 
آلودگی در یک بعد حل می‌شود و سپس نتایج به‌دست‌آمده 
ورودی مرحله دوم محاسبات (حل معکوس معادله انتقال 
آلودگی) را تأمين می‌کند و نتایج نهایی (تابع شدت آلاینده) 
شاه رم توت 

با توجه به شکل (۱) مدل‌سازی انتقال و انتشار آلاینده در 
رودخانه با فروض تقسیم رودخانه به 0 بازه انجام می‌گیرد. طول 
ی قو وووش ان مد فآ سا شش خر مان و 
که در شکل (۱) مشاهده می‌شود غلظت از بالادست به سمت 
پایین‌دست (غلظت از مکان 20 به م6 حرکت می‌کند و 
ترتیب زمانی در حل مستقیم از )به ا و در حل معکوس از دا 
به :) می‌باشد. در روش مستقیم با داشتن الگوی غلظت رودخانه 
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در بالادست (سری زمانی) و یا الگوی غلظت در زمان ابتدایی و 
در تمام مکان‌ها (پروفیل مکانی) غلظت رودخانه در پایین‌دست 
محاسبه می‌شود. در این روش مدل به‌صورت برگشتی عمل می- 
کند یعنی با داشتن الگوی غلظت در زمان انتهایی (پروفیل مکانی 
در زمان 1). غلظت در زمان ابتدایی محاسبه می‌شود. جهت 
استخراج پروفیل مکانی در زمان انتهایی نیاز به داده‌برداری در 
طول رودخانه می‌باشد که در اين پژوهش به دلیل عدم دستیابی 
به داده‌های واقعی, از حل مستقیم معادله انتقال آلودگی استفاده 


۷ ۱ 
1-1 
ی 3 « 
و 


ماحا 
ص 


ما1 1-0 
بازه مورد نظر رودخانه 


شکل ۱ نیمرخ مکانی غلظت آلاینده در رودخانه در بازه ۰ تا ور 


پایین‌دست را محاسبه کرد. دقت محاسبات این دو روش بستگی 
به شرایط مختلفی دارد که نوع الگو یکی از مهمترین عوامل 
یک‌بعدی به‌صورت زیر بررسی می‌گردد [2]: 


۱( ۸6 +۸0 ت ( و ۰ 
۲ ۵ ۵2 0 

در معادله (۸۰6۱ سطح مقطع رودخانه. (1 ضریب پراکندگی 

در جهت طولی؛ >1 ضریب زوال ماده الا تاه و 5 ترم منبع 
می‌باشد. از آنجا که معادله فوق از دسته معادلات دیفرانسیل جزئی 
بدخیم می‌باشد برای تبدیل آن به دستگاه معادلات دیفرانسیل 
معمولی از روش عددی حطوط استفاده می‌شود که پس از 
گسسته‌سازی معادله فوق به روش تفاضل مرکزی ( 060121 
6 معادله زیر حاصل می‌شود که با پارامتر ,1 نمایش 


داده شده است. 
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1 ۹ 2۸ 


اه اه ك 
درل( 02 بر( 0 ۸۳۹ 


6 
(۳) 


گفتنی است که در این مسئله از ضریب زوال واکنش و ترم 
منبع صرف‌نظر شده است. شرایط مرزی و انتهایی برای حل 
عددی معکوس معاد له جابه جایی-پراک‌ندگی در رود خانه 


به‌صورت زیر در نظر گرفته می‌شود: 


00.0 -0 0-3 

0 

)۳( 0 واشت م ۳ 
62۱ 


00:0 0,00 061 


که (0,) شرط انتهایی يا توزیع مکانی نوده م شاهده شده در 

در تحقیق موردنظر ضریب پراکندگی به‌صورت متغیر 
مدنظر است. در این تحقیق از مدل‌های محاسبات جریان که در 
آنها از معادلات سنت ونانت استفاده شده است برای محا سبه 
سرعت. عمق. مساحت جریان و غیره با ضرایب متغیر استفاده 
می شود. روش‌های مختلفی برای تعبین ضریب پراکندگی طولی 
وجود دارد. از جمله این روش‌ها می‌توان به استفاده از روابط 
تجربی ارائه شده توسط محققان مختلف و انجام آزمایش میدانی 
توسط ماده ردیاب و روابط مختلف دیگر اشاره نمود. 

پس از به دست آوردن پارامترهای جریان» به حل مستقیم 
معادله جابه‌جایی -پراکندگی به‌وسیله روش عددی خحطوط ( 1۲6 
دمص گم 0مطاعهه [هع۲0۳۵) و طرح حافظ گروه پرداخته می - 
شود. روش عددی خطوط مفهوم ساده‌ای دارد که در واقع دستگاه 
معادلات دیفرانسیل جزئی را به غیر از یک متغیر مستقل گسسته- 
سازی می‌کند. تون نیمه کشسببته‌سارم دستگاهی از معادلات 
دیفرانسیل معمولی را نتیجه می‌دهد که به‌صورت عددی می‌توان 
آنها را حل کرد. پس از جداسازی معادله (۱) که با روش نیمه 
تسار انجام می‌شود. حاصل مها زا ۱ کست کاه 


معادلات دیفرانسیل معمولی می‌شود که می‌توان با به‌کارگیری 
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طرح حافظ گروه به حل این معادله پرداخت. طرح حافظ گروه 
از فضای افزوده به نام فضای مینکوفسکی (۷110107/5161) استفاده 
می‌کند. تفاوت اساسی مابین طرح حافظ گروه و روش‌های عددی 
دیگر این است که در روش‌های قبلی از فضای اقلیدسی "1 
استفاده می‌شود در حالی که در طرح حافظ گروه از فضای 
مینکوفسکی "۷1۳ استفاده می‌شود. در واقع بعد دستگاه در فضای 
میتکوفسکی یک مرتبه افزانش بیدا مین کند, یکی از مرآبای ابتفادة 
از فضای مینکوفسکی در طرح حافظ گروه اجتناب از جواب‌های 
جعلی (وامتتستام6) و تقاط ثابت غیرواقعی رعاصتمم 1260 ]عمط0)) 
است. 

طرح گروه حافظ در واقع گروه تقارنی دستگاه مورد بحث 
را حفظ می‌کند. برای دستگاه معادلات دیفرانسیل « بعدی 


به‌صورت 


۰1۳ 621۳ 0 ی 


۹3 


که در اين معادله ) بردار 1 بعدی است.] یک متغیر زمانی است 
۶ تابعی با مقدار بردار از ن) و.] که میدان برداری نامیده می‌شود. 
برای سیستم‌های معادلات دیفرانسیل غیرخطی, معادلات را در 
دستگاه‌های دینامیکی تقویت‌شده جای داده می‌شود. می‌توان 
معادله بالا را به یک دستگاه دینامیکی افزوده 1+ بعدی تبدیل 
کر د[27] 


1)),1( 


77" 0 
اصا معع اما 
٩:۱0 ۱۳۳) ۱‏ 


0 
۴, 


به‌وضوح اولین معادله در (۵) همان معادله (ع) اما دومین 
تاه اه شب شا وا باس ها افیه. نی اا ۳ 


2 22 20 (0 


صدق می کند که 2 یک متر مینکو فسکی ( ۷۱0۱۷۹10( عتتاعو) 


است و به صورت زیر نشان داده می‌شود [27]. 


1 1۹ 
8 ۹9 


نشريهةُ مهندسی عمران فردوسی 


حل معکوس (در زمان) معادله انتقال آلودگی در رودنحانه... 


در معادلات فوق ,1 تابع همانی از مرتبه 9 و 1 نمایانگر 
رازه ی 

معاذله (۱) داسته جواب در فضسای منک فسکی را تین 
می‌کند که نت جه آن مخروط پوچ (معمع [[۲9) ش کل ۲ 
می‌شسود [27]. با در نظر گرفتن زمان مثبت در حل مسستقیم و 
زمان منفی در حل معکوس مخروط به دو سمت مثبت و منفی 
که ام طلاحاً به آنها مخروط آینده (606 ۳000۲6) و مخروط 
گذشته (6086 ۳2۵۲) گفته می‌شود تبدیل می‌شود. شکل (۲) 
شرط مخروط را نشان می‌دهد که تاکنون شرط مخروط 
بدیهی‌ترین محدودیتی است که بر روی دستگاه دینامیکی اعمال 
شده است در واقع می‌توان گفت این شکل محدوده جواب 
غلظت,رودشانه در قضای متحوفستگی واشان می‌ده: 


۱ 


ااگزاه 


شکل ۲ شرط مخروط و فضای مینکوفسکی افزوده [19] 


با تعریف ماتریس ۸ به صورت زیر 
20,0 
| ۳" 


۸ :-- ۸ 
21), 0 


گر 


و ثابت در نظر گرفتن آن مقدار آن در گام ها و ی 
با نوشتن نگاشت نمایی برای ماتریس ()۸ در فاصله 
,] +]>] > را هنگامی که ,] <] است رابطه زیر حاصل می‌شود 
[27: 
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1[ ۵ 1 
۰ 8 ی 
9 5 [(۸ اوه 
.2 9 
1 ۳ 


که با استفاده از رابطه (رایر))1< 1 و حایگذاری [([۵ ۸۵| 
به‌جای ([)0 در رابطه 00(0 < بت که یکی از فواید طرح 
حافظ گروه است. معادله زیر حاصل می‌شود [27]: 
|6۱ +6 ۶( ره) 
(۱۰ 1 ۳ ف از 
| ۶ 
در م عاد له فوق ) غاظت رود خانه در ز مان ] و ,0 
غلظت در زمان 1+ می‌با شد. ,9,۰۵ نیز به صورت زير تعریف 
می‌شوند و 1 با استفاده از روش عددی خحطوط محاسبه می شود 
که نتیجه آن در معادله (۲) نمایش داده شده است[27]. 
|| ۱۶| ۸۱ 
ا| 6 || 
ا| ,۵۱۱۴ ر. 
)طصله ع: ط 
اا ب || 
برای حل معادلات دیفرانسیل معمولی به‌صورت معکوس با 
استفاده از طرح معکوس حافظ گروه. معاد له 63 با گام های 
دانت‌گاه معادلات دیفرانسیل ظ بعدی به صورت معکوس معادله 


(4) به شکل زیر نوشته می‌شود؛ 


) 0۱۱) 


2 0۹9 


) 0۱۲ 


ات 


]۲ *۰۰ع ) 


۱۳ 


با یک مقدار نهایی ) مشخص در 0< 1. هدف بازیابی 
مقادیر گذشته ) در زمان گذشته 0 >1 می‌باشد. همانند روش 
مستقیم. معادلات در سیستم‌های دینامیکی تقویت‌شده جای‌داده 
می‌شود. می‌توان معادله بالا را به یک دستگاه دینامیکی افزوده 
1 بعدی تبدیل کرد 1271 


(ع1 


0 
٩ 6 0 | 6‏ 
(۱ ۳ ۳ د 
از و سع | لاه 


گ 
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۳ 


مانند روش مستقيم ادامه مسیر با گام زمانی منفی برای 
روش معکوس نوشسته می‌شسود که همان نتیجه قبل را می‌دهد 
بنابراین مانند روش مستقیم با این کار می‌توان د ستگاه دینامیکی 
اصلی 9 بعدی در "8 را به دستگاه دینامیکی 01 بعدی در 
۳ متناظر کرد که درنهایت با بسط دادن ماتریس (()۸ و 
نوشتن نگاشت نمایی برای آن. معادله زیر حاصل می‌شود که به 
کمک آن می‌توان به حل معادلات دیفرااسیل معمولی به صورت 
معکوس پرداخت. 


هی( وا 
0[ 
0 ام | ۶ 


در معاد له فوق ,) غاظت رود خانه در زمان) و ,6 
مانند روش مستقیم است. 
الگوریتم کلی حل معادله انتقال آلودگی به‌صورت مستقیم و 


معکوس به‌صورت زیر می‌باشد. 


نتایج و بحث 


سه مثال جهت ارزیابی روش مستقیم و معکوس حافظ گروه 
طراخی شده است. این سه مثال به‌خوبی دقت این روش را نشان 
می‌دهند. با توجه به اينکه در این پژوهش معادله انتقال آلودگی با 
ضرایب متغیر حل شده است. ابتدا باید پارامترهای جریان محاسبه 
شود که این پارامترها با حل معادلات جریان و با استفاده از مدل- 
های موجود به دست آمده است. مشخصات رودخانه با توحه به 
مثال‌های واقعی در نظر گرفته شده است و با استفاده از مدل‌های 
محاسبات جریان پارامترهای دبی. عمق» سرعت. مساحت جریان 
و دیگر پارامترها در نقاط مختلف رودخانه به دست آمده است. 
روش‌های مختلفی نیز جهت تعیین ضریب پراکندگی وجود دارد 
که در این پژوهش از روابط تجربی موجود استفاده شده است. 
جدول (۱) مشخصات و سایر پارامترهای هر مثال و شکل (4)؛ 
(0) و (1) سرعت. ضریب پراکندگی و مساحت جریان نسبت به 
طول رودخانه را نشان می‌دهد. 


نشریه مهندسی عمران فردوسی 


1 حل معکوس (در زمان) معادله انتقال آلودگی در رودحانه... 


ضریب پراکندگی در نقاط مختلف مساحت خیس شده در نقاط مختلف سرعت جریان در نقاط مختلف 


استخراج معادله انتقال آلودگی با ضرایب متغیر در رودخانه 


گسسته‌سازی معادله فوق با استفاده از روش عددی خطوط و تبدیل آن به دستگاه 
معادلات دیراد نسی معمولی 


تشکیل دستگاه دینامیکی افزوده و حل آن با گام زمانی مثبت و منفی حهت بدست 


آوردن معادله کلی 


تعیین شرط مخروط و دامنه جواب معادله با حل ماتریس 
0 -< ءج 37ع 


استخراج معادله کلی برای حل معادله انتقال آلودگی به‌صورت مستقیم و معکوس در 
زمان در رودخانه با استفاده از طرح حافظ گروه 


شکل ۳ فرایند مدل‌سازی عددی انتقال آلاینده در رودخانه با ضرایب متغیر 


جدول ۱ مشخصات و پارامترهای هر مثال 


پارامتر / 
و طول دبی الگو آلاینده ضرایب هدف 
مثال | لا | تهج 
مثلث دو قله م‌ تقیم / 
مثال اول ۳۰ ۱۵ متیر ِ حل 
(نامنظم) بازسازی شدت الاینده و پروفیل مکانی در پایین‌دست 
تمانی ده قلة و 
مثال دوم ۳۲۰ ۱6 متیر ۲ 
(نامنظم) بازسازی ستالنت: الاینده و پروفیل مکانی بو بالادست 
مثلث دو قله کت 
مثال سوم ۳۲۰ ۱6 متغیر ۳ 
(نامنظم) بازسازی شلانت الاینده و پروفیل مکانی در بالادست 
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دس( 


۷/۷ 
6 و و و و و وه وه و 
بر سر ها مه نب چم ما چا دا دز تاج 


1 0.80.9 0.60.7 0.40.5 0.3 0.2 0.1 0 
1 
شکل ۶ تغییرات سرعت در راستای طولی رودخانه 


بب 


حرط 
ود اه اه اه هه ها اه 
سم سر ۱ ۵0 لد ص لا ج دیا با تسا 8 


1 0.80.9 0.7 0.6 0.40.5 0.3 0.2 0.1 0 
ب/ 


ابتدا با حل مستفیم طرح حافظ گروه با استفاده از سری 
زمانی به‌منظور شبیه‌سازی روش مستقیم و نشان دادن دقت عمل 
این روش از شرایط مرزی غلظت ورودی به رودخانه مطابق 


هو 
سح ال نظر 
کت 


شکل )۷ و شرط مرزی غلظت خروجی 0 - 
گرفته شده است و پس از شبیه‌سازی برای ارزیابی این روش از 
حل عددی تمام ضمنی (1۳00116010 ولد تفاضل محدود ( ۳106 


6 استفاده شده است. 


جهت حل معکوس معادله انتقال آلودگی نیاز به پروفیل 
طولی یا مکانی است که برای به دست آوردن این پروفیل باید 
غلظت آلاینده در بازه‌های مشخص در یک زمان خاص از 
رودخانه به دست آید که در این پژوهش به دلیل عدم دستیابی به 
داده‌های واقعی, با استفاده از حل مستقیم پروفیل مکانی در زمان 
انتهایی به دست می‌آید. در واقع خروجی مدل مستقیم به‌عنوان 
ورودی مدل معکوس در نظر گرفته می‌شود و به‌منظور صحت- 
سنجی مدل معکوس. پس از شبیه‌سازی مدل معکوس و بازسازی 
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3 


ب 


«جج(]/2] 
ار ی ان نج اب دس نز 
سم سس ص 00 بل ص ‏ خد دنا با سا و 


1 0.80.9 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0 
1 
شکل ۵ تغییرات ضریب پراکندگی در راستای طولی رودخانه 


ض( 


هصت 0 
+ هو و و و و و ه هو و 
سل سر ها مق ید چ ما چا دا دز سا 


1 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0 
هصا/1 


شکل ۷ منحنی سری زمانی غلظت آلاینده در مرز ورودی 


پروفیل مکانی در زمان‌های مختلف. نتایج حاصل از مدل معکوس 
با مدل مستقیم مقایسه می‌شود. با توجه به اینکه برداشت هرگونه 
داده‌ای بدون خطا امکان‌پذیر نیست برای محاسبات و اندازه‌گیری 
شدت آلاینده در این مثال‌ها میزان حطاهای مختلفی در نظر گرفته 
شده است. این خطاها به نتایج حاصل از روش مستقیم که ورودی 
روش معکوس است اعمال شده تا پاسخ‌های به‌دست‌آمده نزدیک 
به واقعیت باشند. به همین منظور خطای ۵ تا ۱۵ درصد به پروفیل 
مکانی غلظت آلاینده در زمان انتهایی اعمال شده است. در انتها 
برای مقایسه این دو حل از شاحص‌های خطای آماری نظیر 
ضریب همبستگی (87) جذر میانگین مربع خطاها (1315۳۲) 
میانگین خطای مطلق (۷1۸2) و درصد خطای نسبی (۷15۲ 


استفاده شده است. 
مثال اول: کاربرد طرح مستقیم حافظ گروه در رودخانه 1 


حافظ گروه طراحی نله استت: در این مثال رودخانه پس از 


نشریه مهندسی عمران فردوسی 


3 


تقسیم به ۳۹ بازه ۰۰ متری و اعمال غلظت ورودی در مرز 
رودخانه, غلظت در بازه‌های مختلف محاسبه شده است. در این 
مثال مدل در چهار زمان بی‌بعد شده (۰,۱۰ (۰,۳۳ (۰,۸۳) و 
(۱) شبیه‌سازی شده و نتایج در شکل (۸) نمایش داده شده است. 
همچنین نتایج ارزیابی با حل عددی تمام ضمنی تفاضل محدود 
با استفاده از شاحص‌های آماری در جدول (۲) آورده شده است. 

با توجه به شکل (۸) دقت طرح مستقیم حافظ گروه در حل 
معادله انتقال آلودگی با ضرایب متغیر در رودخانه بسیار خوب 
بوده و همان‌طور که در شکل (۸) مشاهده می‌شود این حل انطباق 
نسبتاً کاملی با حل تمام ضمنی تفاضل محدود دارد. قله‌ها به خوبی 
بازسازی شده‌اند و باگذشت زمان نیز خللی در بازسازی قله‌ها در 
پایین‌دست رودخانه ایجاد نشده است و انطباق این دو روش در 
تمام زمان‌های شبیه‌سازی حفظ شده است. همان طور که در 


1.1 
۲۱0۱/۵6 ۱ 
۳ 09 
08 
1/۲:::3 (ظ 0.7 
06 
05 
04 
03 
02 
0.1 
0 
1 0.9 0.70.8 0.6 0.40.5 0.3 0.2 0.1 0 
ب/ 
اس[ 
1۲۱۱0۱۶ ۱ 
۶۳5 09 
0 
آحدده)/۱ )4 0.7 
06 
05 
04 
03 
02 
01 


0 
1 0.80.9 0.7 0.6 0.40.5 0.3 0.2 0.1 0 
ن/ 


۰/۵ ۱۱:۹ 


م0 


حل معکوس (در زمان) معادله انتقال آلودگی در رودحانه... 


جدول (۲) مشاهده می‌شود شاخحص آماری نیز دقت بالای این 
روش را نشان می‌دهد. مقادیر بسیار کم 1۸ 451 و ۱1 
و همچنین میل کردن مقدار 18 به ۱ نشان‌دهنده دقت بالای طرح 
مستقیم حافظ گروه و تطابق نسبت کامل این روش با روش تمام 
ضمنی تفاضل محدود می‌باشد 


جدول ۲ شاخص‌های آماری مثال اول 
شاخص آماری مقدار 
)۱۵۸۲ 0.05 
12690 99997 
)۳۱۷۲۹۳۴ 002 
(1۷]۳)6 023 
1.1 
6 ۲۱۱0۱ ۱۱ 
۳ 09 
05 
21-::/) (2 0.7 م 
06 2 
05 6 
4 3 
03 
02 
01 
0 
1 0.9 0.8 0.60.7 0.5 0,4 0.3 0.2 0.1 0 
5/17 
اسب 
۲۱۱۱۸۱۲۶ | ل 1 
0۳ 09 
0 
۷/۸3 0 0.7 ب 
2 
كِ 
0 


1 0.80.9 0.60.7 0.40.5 0.20.3 0.1 0 
ب/ 


شکل ۸ مقایسه طرح مستقیم حافظ گروه با روش تمام < صمن تفاضل محدود 
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مثال دوم: کاربرد طرح معکوس حافظ گروه در رودخانه با 
اماهه از پروفل مکالی:. با ترجه به اینکه ثبت داده‌ها هخشه 
دارای خطا است. بنابراین بر روی داده‌های ورودی مدل معکوس 
(پروفیل مکانی) ۵ تا ۱۵ درصد خطا به‌صورت کاملاً تصادفی 
اعمال شده است تا بتوان اثر این خطاها را در پاسخ مشاهده کرد. 
در این مثال پروفیل مکانی ورودی مدل معکوس با استفاده از 
شرایط ابتدایی (لهثانع؟ مهمتانهه0ع) در شکل )٩(‏ به‌دست آمده 


است و شرایط مرزی ورودی 20 (060,۷ و شرایط مرزی در 


0 ۳ 
در نظر گرفته شده است. در این مثال دوم نتایج 


4 


انتها و - 


با 


در سه زمان بی‌بعد شده (۲۵, ۰ (۰,1) و (۱) اه مت استنت ( 
نتایج بدون در نظر گرفتن خطا در شکل (۱۰) و پس از اعمال 
درصد خطا ۸۱۵ بر روی ورودی مدل در شکل (۱۱) نمایش داده 
شده است. مدل با استفاده از شاخص‌های آماری با حل مستقیم 


ارزیابی شده که نتایج در حدول (۳ ارائه شده استا: 


۳" 


1۱۳۱13۳۰۰ هه 
۳۶۲۶ 


«( ۳۰۰-025 


0/0 
5 9996و و 
سدت ها و6 ید ماج دا دا تام 


1 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0 
ن[/ 


۳ 


010 
اس سس سا سس ۳ 
سس وا مق ید با چد دا دا سا 


۱ ۵ تست ست مت 
۳۶ 


۱( ۰۵۸۵6 


1 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0 
ن/ 


تب 


۱ 
تاکز 


امس و 


هس 
ماه هام اه ماه 
باس ۱۵ 00 لد ۵ را ج دیا ۱۵ سا ق 


1 0.9 0.8 0.7 0.6 0.5 0.4 0.3 0.2 0.1 0 
ن1/< 


شکل ۱۰ مقایسه طرح معکوس حافظ گروه با روش مستقیم 
بدون اعمال خطا 


سال سی و پنجم» شمارة چهار» [ ۱۶۰ 


۷ 


جدول ۳ شاخص‌های آماری در حالت‌های مختلف مثال دوم 


خطا/ شاحص 


((۵/ع)ظ۱۷۲۸ 
(12)00 
(0۵/۲۲ 1۱۷/5۲ 
107 


601001 
0۹6 


019 
9962 
008 
234 


601001 
570 


038 
97.1 
017 
639 


6۲۲0۲ | ۲ 
157 10۳۹7 
017 039 

94.04 ۱ 9 
039 025 

12.56 | 5 


قت 


مان 
+ هه مه هي ه ه ه ه ۰ 
با صا 00 لد ص با خجد تا با تسا 


0 0.1 0.20.3 0.40.5 0.60.7 0.80.9 1 


ن 1/1 


شکل ٩‏ پروفیل مکانی (شرایط ابتدایی) مدل مستقیم 


۱-۰ 1۱ بت 2 ۳ 
5 13۰ 


2( 1/۰۲05 
۰۳۲0۱154 


ن/*< 


۱/۸۱6 (ظ 
۰۲۲0۱-215 


0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 


(‌ 


0/0 
99 و و 
سس ها م6 ید ما چا با دز سا 


۳ 


مسا 
:+ 99ووووووو 
سس ها م6 ید ه ماج با دما 


0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 


نآ 


1۱: 


ات گر 


دسا ۱ ۵ 
56 ۰۲۳۲0۱ 


۳" 


م0 
660 9966و 
ساب ۱۵ 00 لد ض زا جح دنا ۵ سا 8 


0 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1 


نآ/< 


شکل ۱۱ مقایسه طرح معکوس حافظ گروه با روش مستقیم 


با اعمال خطا ۱۵/ 


نشریه مهندسی عمران فردوسی 


1۸ 


در اين مثال چون مدل مستقیم با استفاده از شرط ابتدایی 
اجرا شده است مدل معکوس محدودیت زمانی در اجرا مدل 
ندارد و بازسازی مدل بازگشتی به‌طور کامل صورت گرفته است. 
همان طور که از شکل‌های (۱۰ و ۱۱) پیداست هر چه زمان 
بیشتری از شبیه‌سازی می‌گذرد دقت محاسبات کمتر می‌شود با 
ان سعال همان طون که ششاهله فی شود اقله‌ها با دقت سیای غویی 
بازسازی شده‌اند. دقت مدل می‌تواند به عواملی همچون تعداد 
گام‌های مکانی و زمانی. زمان شبیه‌سازی, از بین رفتن اثر الگوی 
ماده آلاینده در رودخانه بستگی دارد که به دلیل تعداد کم گام‌های 
زمانی و از بین رفتن اثر دو قله غلظت با گذشت زمان نتیجه بسیار 
خوبی حاصل شده است. یکی از نقاط قوت این مدل عدم نیاز به 
تقسیم‌بندی رودخانه به تعداد بازه‌های زیاد می‌باشد که به دلیل 
مشکل بودن برداشت داده در طول رودخانه این یک مزیت مهم 
محسوب می‌شود. همان طور که از مقایسه شکل‌های (۱۰ و ۱۱) 
مشاهده می‌شود با افزایش خطای اعمال‌شده بر داده‌های ورودی 
مدل معکوس میزان دقت بازیابی الگوی آلاینده کم می‌شود. 
مقایسه سنون‌های جدول (۳) نیز نشان می‌دهد با افزایشس خطا 
مقادیر 0۷۸ 55 و 118۳ افزایش پیدا می‌کند و از میزان 
8 نیز کم می‌شود که این تغییرات نشان‌دهنده کاهش دقت 
بخاشیات. با افرایشن. طا می‌باشد. همین با اغمالن شطا: 
نوسانات غلظت در پایین‌دست رودخانه در زمان‌های انتهاپی 
شبیه‌سازی مشاهده می‌شود که با افزايش خطا این نوسانات بیشتر 
می‌شود. برای نشان دادن دقت روش نتایج مقایسه با شاحص‌های 
آماری در جدول (۳) آورده شده است. همان‌طور که در جدول 
(۳) مشاهده می‌شود شاخحص آماری 1۷۲۵۳ 2۷۲۹۲ و ۷1۲ 
تقاانی کی نلک هف وی یو شقن دا فان ما 
همچنین مقدار بالای **1 نشان‌دهنده تطابق خوب این دو روش 
می‌باشد و حتی این مقادیر هنگام اعمال خطا نیز نتایج خوبی را 
نشان می‌دهد که بیانگر این است که در صورت به وجود آمدن 
خطاء الگوی آلاینده به‌حوبی بازسازی می‌شود و نتایج بسیار قابل 


قبولی به دست می‌آید. 


مثال سوم کاربرد طرح معکوس حافظ گروه در رودخانه با 
استفاده از سری زمانی. با توجه به الگوی سری زمانی غلظت 
آلاینده (غلظت در مکان ابتدایی رودخانه و در تمام زمان‌ها) شکا 


(1) و شبیه‌سازی مدل مستقیم. پروفیل مکانی در زمان انتهایی که 


نشريهةُ مهندسی عمران فردوسی 


حل معکوس (در زمان) معادله انتقال آلودگی در رودنحانه... 


در واقع ورودی مدل معکوس است به دست می‌آید. در اين مثال 
نتایج در سه زمان بی‌بعد شده (۰,۳۵ (۰,77) و (۱) آورده شده 
است. در اين مثال رودخانه مورد نظر پس از تقسیم به ۳۹ بازه 
۰ متری و اعمال پارامترهای جریان و ضرایب متغیر» شبیه- 
سازی عددی مدل انجام گرفته است و پس از آن نتایج مدل که 
با درصد خطاهای مختلف به‌دست‌آمده در شکل (۱۲) و (۱۳) و 
سپس مدل با استفاده از شاخص‌های آماری با حل مستقیم ارزیابی 


شده که نتایج در حدول 613 ارائه شده آمخستا: 


جدول ۶ شاخحص‌های آماری در حالت‌های مختلف مثال سوم 


۰101 ۰01 ۰01 ۰01 


خطا/ شاخحص 
ِ 12 59 106 130 


0.91 073 056 0۳26 ۱۲۷۲۸۲ 


92.09 | 9481 | 95.14 2 12)90( 
05 031 024 0.11 ۱ ۳۱۹۲) 
113 102 52 46 (۷۲۳۳)۶۵( 


علی رغم اينکه زمان مدل مستقیم بیشتر از مدل معکوس 
است. مدل معکوس تنها توانایی شبیه‌سازی نیمی از زمان بوده 
است یعنی آن قسمت از رودخانه که غلظت از مرز ورودی وارد 
رودخانه نمی‌شود بازسازی برگشتی صورت می‌گیرد که در واقع 
مهم‌ترین محدودیت این مدل هست. همان طور که از شکل (۱۲ 
و ۱۳) پیداست روند بازگشتی به‌حوبی صورت گرفته است و 
قله‌ها با دقت خوبی بازسازی شده‌اند. همان طور که از شکل‌های 
(۱۲ و ۱۳) پیداست هر چه زمان بیشتری از شبیه‌سازی می گذرد 
دقت محاسبات کمتر می‌شود که این کاهش دقت هنگام تشکیل 
قله‌ها در زمان انتهایی نیز معلوم است. با توجه به شکل (۱۲ و۱۳) 
کاملاً مشهود است که "هرچه درصد خطا اعمال‌شده بر روی 
پروفیل مکانی بیشتر شود بازسازی مدل برگشتی با دقت کمتری 
همراه خواهد شد. کاهش مقدار ۸۲ 31558 و 158 و 
افتایشن تقدار نوی ات اس رشان می‌دهد. گفتتی اشتت 
این مقادیر هنگام اعمال خطا نیز نتایج خوبی را نشان می‌دهد که 
بیانگر این است که در صورت به وجود آمدن خطاء الگوی آلاینده 
به‌حوبی بازسازی می‌شود و نتایج بسیار خوبی به دست می‌آید. 
همچون مثال قبل با اعمال خطاء نوسانات غلظت در پایین‌دست 
رودخانه در زمان‌های انتهایی شبیه‌سازی مشاهده می‌شود که با 


افزایش حطا این نوسانات بیشتر می‌شود. برای نشان دادن دقت 


سال سی و پنجم» شمارة چهار ۱۶۰۱ 


امیر محمد سعادت- مهدی مظاهری - جمال محمدولی سامابی 


روش نتایج مقایسه با شاخص‌های آماری در جدول (4) آورده 
لو اتجنت: همان‌طور که در حدول (۶) مشاهده می‌شود شاخحعص 
آماری 1۷۲۵۳ :18151 و 1158 مقادیر کمی دارند که دقت خوب 


روش معکوس را نشان می‌دهد همچنین مقدار بالای 12 
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شکل ۱۳ مقایسه طرح معکوس حافظ گروه با روش مستقیم 
با اعمال خطا 7/۱۵ 


سال سی و پنجم» شمارة چهار» [ ۱۶۰ 


1۹ 


نشان‌دهنده تطابق خوبت این دو روش می‌باشد. افزایش ص1۷ 


5 و 18۳8 و کاهش 8 براثر افزایش خطاء کاهش دقت 


مدل براثر افزايش خطا اعمال‌شده بر شرط ابتدایی مدل معکوس 


۱۳۰۰( هس بت ست سا 


۳۶:۳۹ 


2( 0:۸۸:25 
6۲۲0۱-۵ 


0 0.1 0.2 0.3 0.40.5 0.6 0.7 0.80.9 1 


1۱(۲ دس بت بت سر 


۳۶۹ 


0( 0/۸: 
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1 0.80.9 0.7 0.6 0.40.5 0.3 0.2 0.1 0 
ن/ 
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۳۶۲۹۵ 
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0۲۲۵۲۰۵ 
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هس01 
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فصست0 
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9/0 
و همه هه ۵ :2:6 
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شکل ۱۲ مقایسه طرح معکوس حافظ گروه با روش مستقیم 


نشریه مهندسی عمران فردوسی 


نتیجه گیری 
در این پژوهش به‌منظور به دست آمدن شدت آلاینده از حل 
عددی معادله انتقال آلودگی با ضرایب متغیر در رودخانه استفاده 
شد. سه مثال جهت ارزیابی این روش طراحی گردید که در مثال 
اول با استفاده از شرط اولیه (مثلث دو قله) دقت مدل مستقیم 
مورد ارزیابی قرار گرفت که نتایج به دست آمده نشان داد روش 
مستقیم تطابق بسیار زیادی با روش تمام ضمنی تفاضل محدود 
دارد. در مثال دوم به ازای شرط اولیه و با استفاده از روش مستقیم, 
پروفیل مکانی به دست آمد که به‌عنوان ورودی مدل معکوس 
استفاده شد. در این مثال نیز شبیه‌سازی با اعمال درصد خطاهای 
۵ ۱۰ و ۱۵ درصد بر روی پروفیل مکانی صورت گرفت و نتایج 
حاصل نشان داد که قله‌ها به‌طور تقریباً کاملی تشکیل شدند که با 
توجه به از بین رفتن اثر دو قله. نشان از قدرت و توانایی مدل 
معکوس در تشکیل الگوی آلاینده می‌دهد. از خروجی مثال اول 
به‌عنوان ورودی مثال سوم استفاده شد که در این مثال با استفاده 
از طرح معکوس حافظ گروه بازسازی مدل برگشتی انجام گرفت. 
شبیه‌سازی صورت گرفته با اعمال درصد خطاهای ۵ ۱۰ و ۱۵ 


وت 


حل معکوس (در زمان) معادله انتقال آلودگی در رودنحانه... 


بر روی ورودی مدل معکوس انجام شد که در تمام حالات قله‌ها 
با دقت بسیار خوبی تشکیل شدند و نتایج شاحص‌های آماری 
دقت این روش را نشان داد. با توجه به این سه مثال دقت خوب 
طرح حافظ گروه در تشکیل تابع مت ال بنده در وو ان و 
تااشازی توف سکانی فان هساو مد سین 
امکان نشان دادن پروفیل مکانی غلظت در زمان‌های مختلف که 
کمک شایانی در مشاهده تحولات و تغییرات مکانی توزیع غلظت 
ی کید را دار استت ان زرشن سکلانی عون اشتمای و مد 
سازی ترم پایداری که به دلیل عددی بودن این روش است را 
ندارد که همین دلیل باعث بالا بردن سرعت و کم کردن خحطای 
احتمالی می‌شود. علاوه بر موارد فوق یکی از مزیت‌های دیگر 
این روش عدم نیاز به تهیه داده در تعداد بازه‌های کم در طول 
رودخانه می‌باشد که به دلیل دشوار بودن این کار یک مزیت مهم 
محسوب می‌شود. در نهایت به دلیل سرعت و دقت زیاد این مدل 
و ساده بودن محاسبات باعث صرفه‌جویی در زمان می‌شود. از 
محدودیت‌های این روش نیز می‌توان به عدم توانایی در بازسازی 


سری زمانی اشاره کرد. 


۱۱۳۹6۱۱۱۱۱۵۱ 9 


۰ ۷0۲۱ ها ,تمطوتاهانظ ۱۷۵/۵۸۰۲۱11 ۱0۵018۵ ]۲۷۵۵۲-۰۵۵ ۱۵6۵و ب0۵01۲۵) 2۰ ٩۰‏ [2] 


۵ ط10) ]۲۵ ب۶۵ م۷۱۵0 اهمتاممصم‌طه]۷ رتصعصهه ۷۰ ,۱۷۲ ,۲۲ مصه تحص نلج۲۷ 0مصصحصهطه]۱۷ .ر رتتع‌طه۱۸2 ,۷۲ [3] 


۰ ,310-321 .00 ,16 ۷۵۱۰ ,۳0۲۵ ۱۱۷۰۵ روته 1۷ صا 0611۴1621101[ 


۵۲ ۵0۵ ۵0۶ موی ممناهم10020)1 )فصو وا رهتیمصو مه نک بقلاعه))ظ .0.۳ ,0فماوممماجهد رنایه بملتمل [4] 


0 60۱6۲۵۱۱66 1۱016۲۳۱۱۵۱۱۵۱۱۵ 51 - و0 00/۵۲۵۸۵۵۵۰ ظ ی۵۷۵۲ 0109800601 1۳072017۵ صح وه ممتاهصتطصهاطوم 


2)0015(۰) رکا00ع ی 000 و0۱۱ ۱و ۸۹۵۱۱۵۵۱۱ ۲۱۱۷0۱۳۵۱۱۱۵۵ 


0 ۲۱۷۵۲۶ اطونهاو امصصه صا کاصمل تمصا ممتتااآمص مج ۶۵۶ [0مصظ متاهع) ص106 ممسنامو اصهاتاوظ رت ,26 200 عصهدا2 .و [5] 


۰ ,1955-1963 00۰ ,4 .20 ,7 ۷۵۱۰ .56 ۲۷۵/۵۲ ۸۳۸۵۱۰ محططات0ع 21 متاعصهع وم 


۲ 501۷1۲2۵ 10۲ ولمطع]۱۷ نع اهعصاطازتمعاه )0۵۵ ها رلعصهطم۱۷ .و فصج (تقمصییک) 0۰ ۱ [60] 


۰ 20[ ,45-57 .00 وا .20 و131 ,۷۲۵۱ .۷۵۱۱۵8 ۱0 کم ۲۷۵۸۵۲ ۰۷۰ رکصطع۳۲۵9۱ عورعبصا صمتا1060111162 عنام 


ممتاماصسصتی مملصنا 0موج-ه0 مصلونا وممامه ممتاباامظ مهلصامتن گ۵ ممتامه‌گتامم0] رصتطدظ .1 فصع صعطام]۱۷ نک ٩.‏ [7] 


۰ ۷]2۳[ ,101-109 00۰ ,2 00۰ ,11 ۷۵۱۰ ,۵ ۵ ۰۷۰ ,۱۷۵001 طمتاج2تصتاون0 


کصتمج ممتاطامم آه ممتاه‌تتاجم10 مصصتا معوعام۲ قمع مم2عم رتصقصصهی ۷۵1۲ 0محصصصنهطم۱۷ .ل فصع رتتعطه۱۷۲22 ۱۲۰ رقصقطت ,۸ [8] 


164-1714 .00 ,180 ۷۵۱۰ ,۱۵۱۱۵۵6 ۱۱۷۱۲۵۰ ۱۰ و۱۵۱۵ اتازماجا۳۳۵ 0107270ظ مها رم 0عقوه 0۵0۷۵۲۲65 ۲۱۷۵۲ طا عمتنا80 


2016۰. 


م1062 همه چم 0عقهها وتماه ومواه طا ما6۵ کصلممر اصعصتصهاومه گه صمتاه‌گتاجم0] رهز ۷۰ 28 عصفطن ۴۰ ۷۷۰ [9] 


نشريهةُ مهندسی عمران فردوسی 


سال سی و پنجم» شماره چهار ۱۶۰۱ 


امیر محمد سعادت- مهدی مظاهری - جمال محمدولی سامانی 3 


۰ ۸0۲ ,397-410 .0۵ ر4 .20 ,33 ۷۲۵ ۵۵۵ ۲۷۵۸۵ ۸0۲۰ ۰ ب0مطاممود صمتاعصن ماتعصمل اتازماهما10م 

ما ط0امصتصمادمه ۵۴ فصمتاهتعوهاه مامت)باهه عصتعدا معط وهتازجاهمامر 020 موظ رجمع۱۷۷1۱ با .ل 24 تمننچم‌نام۳۲ ,۱ :16 [10] 
۰ .۵ ۵50۱۲ ۲۷۵۸۵ رط۹۵96۲۷۵10 وماتان۱ واامونطمعوعه۱ مطا که وممنامو ممتافصتصهاوم تماهصامتع بتاصم10 
۰ ۳۵0۰ ,1-14 .0 ,2 .0ظ ,41 

۰ ,۵005 1062010 ممتیامو اصمنطامم ۵۶ ومع ت۷۵ 000 صهو ورماجههاه1؟. رتصتصع2 بخ 0صج یمه 1 .6 ۷۲۰ مقاهمینن ۳۰ [11] 
.۰ ,169-182 .00 ,1 .10 ,29 ,۷۵۱ .کوک او .۵5 2۲۱۷۱۲۵0۲۰ 

02260 و1۳۷۵ ما ممتایتتاصم10 ممتامو ممتتاامن صتمص تمه مهمتجره ۵۷ رعط2 ۷۷ ۲۰ 200 رلصا ,2 و220 .ل ,۷۷202 .1 [12] 
۰ 000 ,739-774 00۰ ,241 ,۷۵ ,۵/۱ ۱۱۷۵۵۴۸۰ و.مطامجط واتاتماهدام تج موه معط 

ممتمنوع ممافممرولل-ه60010 207 معط ۶ ممتامو موعم1۳. رتصمره تله ۷ هام۷ .ز قصه بتتع‌طه۱۷۲۸2 ,۱۷۲ ,1تاهطفتام1 ,۱۷۲ [13] 
(صهتویعظ 12) .2019 ,103-113 .00 با .۱0 ر33.3 ,۷۵۱ .۱8 ۸۸۵6۱ .9۱۵۲۲۱ صماه011)ع10 ععسیامو ممتتامم 10۲ و۲1۷ 1 

عطا ما ممتانامی لحمت وله مصصتاعک1۳۷۵ مج طا ممتتامم ۵۶ مممنامو فص 10۵ رتمطج۱22 ,۱۲ 24 طممصععهیرع۳ ۴۰ [14] 
(صقلک۳۵ 12) .2021 ,25-39 00۰ ,2 .10 6 ,۷۵۱ 0۵/۵( رطمتاهتا60 )0۲حمقصج متام 

-0۱۵51) 0۶ ۱۷۲۵۵۵0 :مصیتاط )صمصتصماممیت ماج لصتام0) و ۵۴ وماعلا۲ مطا عصتنه ۳۴66۵۷ رهاو‌دافک .[ .2 24 ومعهاو ,۲ ۲۰ [15] 
5۰ ,2669-2673 00۰ بأ3 ,۷۲۵۱ عع ۵۲0 ۱۷۵/۵ ,۲۵۷۵51011 

مهو 60ومآمصه صا فصمتاهعم1 ممیامو عصعصتصماومه مامانهم‌تانهم ۵۶ مملامعاصع0]؟ رعصج/۷۷ .5 0ص2 رثن ۲ رعطقط2 .1 [16] 
-667 .00 ,1 ,۷۵۱ , [ 201 ۱۱۵۵ ۳ وخآمت۵ ۲ 0۱۱00۲ ۵۹۵/6۲۵۱۵۵ ۱۱۵۱۵۵۵۵۵ 72۸ 13" رعصتا0۵1مصص مرن 1۳۷۵۲90 
2011(۰) ,672 

6 م13 مصتامنااعم۲۵60 10۲ 05مصامصه موعبصا ۵ ممقانهم‌جطمن. مجموا1 ۷۷ با رل 4ص2 روتعطمت0ظ .ظ وتعناهم‌ناهل ,۷ :3 [17] 
00۰ ,2469-2475 .00 ,9 .10 ,36 ,۷۵ کمک ۵۲0 ۷۵/۵ ما50 صمتامطتطهومه 0۷26۵صامع 2 ۵۴ وتماونط 

,47 ۷۵۱۰ ,1۳۵۹ و۱۷۵ ۵۱ .۱ ۲۱ " رقصجهاهامنج ممتامب0صمی فقمط 0تهعامجه ۶۵۲ مصصمطمو مصامومنم متام رباشناً :5 0۰ [18] 
.4 ,2567-2576 00۰ ,12-13 .10 

فص ما 02۵2 نا نم مصمدامی موم نامع مهن فط ۳1 رعصعطن بظ رل فصع رعصقط ۲۷۷۰ .) پتالناً .9 0 [19] 
۰ ,61-80 00۰ ,1 .10 ,26 ۷۵۱۰ ,5 27۰ م۵۲1۵ ۱/۵۵۵5 ۱۷۸۵۸۸۵۰ رجمتاهناوه حموو 20۷76010108-5 

)-صعک1 تهعمتاصمه لمممتعصمصصتل ما له مهن مصز مهبم مهمتممرج متامتع-علن ۸۵ ومنککا 2 .0 20 عصودن ۷۷۰ 2۰ [20] 
7 ,590-598 .00 ,101۷ ]15 .20 ,79 ,۷۵۱ :8۵ ۳۳۵۵۵0۲6 بطمتاعتاوه 

"م600۵0105 ومهامرق1 معط ۵۶ صماماميم وطمنیهت مطا رم مصمدمه موم متام رتم‌طوه]۲ .۵ ۱۷۲۰ 4ص2 برل‌صح‌دامج‌طاطام .9 [21] 
۰ ,1003-1009 .00 ,6 .10 و35 ,۷۵۱ ,.27 0۱0 ۸۱۵۰ .۲۲۵ 

۴ 000)عظ 20 معصمطمو مصتهومت وتامنع مطا مصتل‌تمعمط ربمت .و 4صه رفاصمافوظ ,۷۲ ۲۰ رتطعنطع۳2:)۵02 ۷۲ رمع ۷۷۰ [22] 
۰ ,102555 .۵ ,15 ,۷۵۱ ,وا کاعع راه0مصه مملتمن مزع ۵ فصمتهآنصتی لهعزرمصصناه مه مع فص 

0ص ونامعصهمعممومطصمه لقممتعصمصتل-ت ناه رم مصمدمی نومیم نامع 0200۳۵۲0 مط ۳1 وتالنا ٩.‏ .0 20 عصعطن ۷۷۰ ۰ [23] 
۰ ,4027-4048 00۰ ,15-16 .۱0 ,38 ,۷۵۱ .۷۵01 ۷۵/۰۰ 41۰ رهم0۲00 ۷۷۵۷۵ 20107/270 تهمصلطمه 

۲ عمصصمطامو ونامزع-ماب ارمنآوه صرح )تمدامرص تملرن قاط تعتقصج ۲ ووه]۱۷ 4صه ۳۱۵۵۲ ۵۶ اقصمل اقصمتاقطصهاص رصهط ۷۰ [24] 
۰ ,1016-1029 00۰ ,101 ۷۵۱۰ ,77۵75۴ کعه۱۷ ۵۵۶ ۱۰ ۱۱۱۰ رقججهاهامتن ممتامن‌جمی اققط 0201022۳0 عصانآ0و 

ونم اطع ما تلهم وتونم مروزنا مطا گم صمتافصصتاو۴؟ راو مصصج۸0 ۲۰ .۲ هه رصهتطلازمطه ۲ رتفمتطت .ظ ۱0۵21 5 [25] 
۰ ,1-12 .00 ,5 .20 ,143 ,۷۵۱ ,2 ۵۲۵۵( ۰ ,۱۱۵061 هنومن تقلنامجن 2 عمزونا 


۶ طموتنومصم 0عوهها-مصناام0مص اهمزمصال. رعصهنعو ۸ 224 رعصه نا .2 ۳۱۵۵520 .۲ وت001 با مل1۵206272 ,۷۲ [26] 


سال سی و پنجم. شمارة چهار» ۱۶۰۱ نشریه مهندسی عمران فردوسی 


۲ حل معکوس (در زمان) معادله انتقال آلودگی در رودحانه... 


-808 .00 ,7 .20 ,64 ,۷۵۱ و.ل .561 0۲۵1 روته ۲۱۷ صا ا0۲حفصقا ماتامو ۶۵۲ فهانترم؟ کصمله‌تللهمه مملوهمونل تمصنل‌بنهم10 
۰ ,919 
,36 ۷۵۱۰ ,۷۵61 ۵۱۱۵۵۰ ۸۷۵۷۸۰ ۲۲۱۰۱۰ رومحجمدمه هطا وفع نامع مج صموری لمعتصحصرل تقمصت-صمه ۵ مصرمن؟ پتاانا ٩۰‏ 02۰ [27] 


20.7, 00. 1047-1068, ۰ 


نشریه مهندسی عمران فردوسی سال سی و پنجم شماره چهار ۱۶۰۱ 


